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Abstract 

This review paper presents the known results on the asymptotics of the 
survival probability and limit theorems conditioned on survival of critical 
and subcritical branching processes in IID random environments. The 
key assumptions of the family of population models in question are: non- 
overlapping generations, independent reproduction of particles within a 
generation, independent reproduction laws between generations. This is a 
biologically important generalization of the time inhomogeneous branch- 
ing processes. The assumption of IID (independent and identically dis- 
tributed) random environments reflects uncertainty in the future (as well 
as historical) reproduction regimes in actual populations. 

This review focusses on a particular range of questions of prime interest 
for the authors. The reader should be aware of the fact that there are 
many very interesting papers covering other issues on branching processes 
in varying and random environments which are not mentioned here. 

1 Introduction 

There has been a constantly growing interest in branching processes in ran- 
dom environment (BPREs) following the pioneering publications by Smith and 
Wilkinsonin g7] in 1969 and Athreya and Karlin [13], [H] in 1971 (see [DD] for a 
list of references up to 1985 and [TJ-[T1], HB]-[3T], [15] -[3D], El]-[37], @3J, [g], 
[48] . and [52]-[64] f° r some more recent results). A BPRE is a stochastic popula- 
tion model where individuals constituting a generation reproduce independently 
according to a common offspring distribution. The random environment feature 
of the model means that the reproduction laws randomly change from one gen- 
eration to the other. In this review we summarize some recent foundings on 
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critical and subcritical discrete time BPREs under a key assumption that the 
reproduction laws of particles of different generations are independent and iden- 
tically distributed (HD). The results reported here are mainly due to the joined 
efforts of the Russian- Germain team of researchers consisting of V. Afanasyev, 
C. Boinghoff, E. Dyakonova, J. Geiger, K. Fleischmann, G. Kersting, V. Vatutin, 
V. Wachtel and others (see also Q2] for an earlier review). 

In the study of stochastic processes in random environment it is crucial to 
recognize the difference between quenched and annealed approaches. Under the 
quenched approach characteristics of a BPRE such as the survival probability at 
time n are treated as random variables, where the source of randomness is due to 
uncertainty in possible realizations of the environment (see, for example, |13] , 
[14], [15], [49], [55] and the relevant literature in [66]). While the annealed 
approach studies the mean values of the mentioned characteristics as a result of 
averaging over possible realizations of the environment (see [5] , [3] , [B] , [H] , [31] ; 
[4"2"] . [35] , [52] . [53"] and the bibliography in [5B]). The annealed approach can 
be viewed as a summary of the more detailed quenched analysis. One has to be 
aware that in biological applications the predictions by the annealed approach 
might be misleading: in reality the environmental history, being uncertain, is 
unique. 

The remainder of the paper is organized as follows. Section [5] introduces an 
important classification of BPREs into three major classes: supercritical, criti- 
cal, and subcritical branching processes. Here we discuss the key assumptions 
imposed on the reproduction laws in the critical and subcritical cases. Section 
[3] contains asymptotic results obtained using the annealed approach. In Sec- 
tion |4] we summarize the results for the multi-type BPREs where individuals 
within the same generation are allowed to have different reproduction laws. The 
quenched results for the critical BPREs are presented in Section [5J The last 
Section 6 is devoted to the reduced critical BPREs. 

Observe that in the setting of random environment we use different fonts for 
the probability and expectation operators depending on whether we condition 
on the environment ("P,£) or not (P, E). Yet another font (P,E) is used when 
we introduce an auxiliary probability measure in the subcritical case (compare 
for example Conditions [5] and [7]) . Constants denoted by c may take different 
values in different formulae (the same reservation holds for c\,C2). 

2 Classification of BPREs 

The basic classification of the branching processes in a constant environment 
recognizes three major reproduction regimes: subcritical, critical and supercrit- 
ical. It is based on the mean value m of the offspring number and these three 
classes are defined by the following relations: m < 1, m = 1, and m > 1. 
Adding the feature of random environment makes the classification issue less 
straightforward due to random fluctuations of the consecutive offspring mean 
values. 
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2.1 BPREs and associated random walks 

A branching process in a varying environment is most conveniently described 
in terms of a sequence of probability generating functions (PGFs) 



i=0 



where qf^ 1 stands for the probability that a particle from generation n — 1 
contributes to the next generation n by producing exactly i offspring (indepen- 
dently of other particles). Time inhomogeneous Markov chain {Z n } describing 
the fluctuations in population sizes is then characterized by the iterations of the 
consecutive reproduction PGFs 

E[^] = .M/ 2 (... (/„(*))...)), (1) 

if we assume that Zq = 1. Taking derivatives in the previous relation gives an 
expression for the expected population size 

E[Z n }=m 1 ...m n , (2) 

where m n = f n {l) stands for the mean offspring number. 
An important example is the linear-fractional case 

f n {s)=r n + (\-r n )- 



1 -(!-*„)« 



fully characterized by a sequence of pairs of parameters (r„, t n ) <G [0, 1) x (0, 1]. 
The iteration of linear-fractional PGFs is again linear-fractional 

MM. . . (/„(,)) ...))= rW + (1 - r in) ) 1 _ {l _ tin ))s 



with 



I _ r (") 

l-t (n) l-t n 1 - 1-ti 

. fn) = —: + —7 m « + • • • + — — m 2 ...TTln. 

I w t n t n -l tl 



mi . . .m n , 



It follows, 



-7-t = 1 + m n + m„_im„ + . . . + mi . . . m„, 

1 - r n 1 - r„_i 



and since V(Z n = 0) = r<"), 

1 



nz n >o) 
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with ak — }_l k and Sk — In mi + . . . + lnm^. The obtained equality indi- 
cates that the asymptotics of the survival probability P(Z n > 0) under certain 
assumptions is governed by the minimal value of consecutive expectations of 
population sizes 

mm(e s \ . • • , e s ") = e ™*(*i, •••,*„)_ 

Notice also that in the linear-fractional case (Z n — 1 | Z n > 0) is geometrically 
distributed Geom(^™^). 

Given a sequence of random PGFs (Iq, i<2, ••-) one can speak of a BPRE 
{Z n}^Lo- Throughout we will assume that 

• the defining random PGFs in are IID so that all F n = F have the same 
marginal distribution, 

• the BPRE starts by a single particle, Zq = 1, unless it is clearly stated 
otherwise. 

For BPREs the relation {1} transforms into 

E | Fi = h,F 2 =h,...]= fx{h{. ■ ■ (Ms)) . . .))■ 

Throughout we will denote the probability and expectation conditioned on the 
environment by 

V(-):=F(.\F 1 ,F 2 ,...), £ [•] :=E[- \F 1 ,F 2 ,...]. 

Observe that the quenched counterpart of the relation ^ can be viewed as a 
realization of a Markov chain £ [Z n ] = e Sn defined by the so called associated 
random walk (ARW): 

S = 0, S n = X 1 + ...+X n , where X n = In F' n (1) . 

Here the jump sizes are IID random variables with 

X n = X:=\nF\l) = \n£[£], 

where £ is a random variable representing the offspring number for a single 
particle reproducing in a random environment. 

The key representation £ [Z n ] = e Sn leads to the following extended classifi- 
cation of branching processes: a BPRE is called (A) supercritical, (B) subcrit- 
ical, (C) non-degenerate critical and (Co) degenerate critical, if its ARW a.s. 
satisfies one of the following conditions 

(A) lim„^oo S n = +oo, 

(B) lim n ^oo S n = -co, 

(C) limsupn^QQ S n = +oo and liminfn^oo S n = — oo, 
(C ) S n = for all n. 

This classification, due to [8], is based on the crucial fact (see for example, 
[5Tj . ch.12, §2) that any IID random walk can be attributed to one of these 
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four classes. For the processes in which E [X] exists this classification coincides 
with the standard classification of branching processes: a BPRE is called su- 
percritical, critical or subcritical if E [X] > 0, = 0, or < 0, respectively (see, for 
instance, [T5J). However, the extended classification seems to be more natural 
and is justified by a series of results published in jS], [30], [55J, [S3], [ST], [55] . 
In particular, in view of the estimate 

V (Z n > 0) = min V(Z k > 0) < min £ [Z k ] = e min <- Sl <-> s ^ 

l<k<n l<k<n 

it follows that the probability of extinction of subcritical and non-degenerate 
critical BPRE equals 1 a.s. Moreover, it was shown in [SJ, [30], [55J, [515], [57] . 
and [55] that the asymptotic properties of the non-degenerate critical BPRE 
with E [X] = are quite similar to those in which E [X] does not exist. 

In what follows we mostly focus on classes (B) and (C) excluding the de- 
generate critical class (Co). Therefore, in the case (C) we will speak of critical 
BPREs omitting the specification 'non-degenerate'. Note that by excluding the 
class (Co) we disregard the critical Galton- Watson processes in constant envi- 
ronments. 

2.2 The Spitzer-Doney condition in the critical case 

We further specify the general condition (C) by assuming that the ARW satisfies 
the classical Spitzer condition with parameter p : 

1 ™ 

- VP(5 fe > 0) -> p e (0,1) as n^f oo. 
n 

fc=l 

According to [23] the Spitzer condition is equivalent to the following Doney 
condition with parameter p 

F(S n > 0) -> p e (0, 1) as n^oo. 

In this text we refer to this condition as the Doney-Spitzer condition expressing 
it in the form 

1 ™ 

p= r lm _Vp (S k > 0) = lim P (S n > 0) € (0, 1). (3) 

n—toc Ji — ' n—too 
fc=l 

It is well-known that any random walk satisfying Q is of the oscillating type 
(see, e.g., Section XII. 7 in [31]). We note that condition ([3]) covers random 
walks satisfying Condition[T](see below) as well as all non-degenerate symmetric 
random walks. In these cases p = 1/2. 

Condition 1 The distribution of X has zero mean and variance < a 2 < oo. 
It is non-lattice. 



5 



Of course, one can consider the lattice case as well. However, this leads to 
unnecessary complications in the statements of the respective results which we 
prefer to avoid. 

Other examples where the Doney-Spitzer condition fulfills are provided by 
random walks in the domain of attraction of some stable law, see Condition [5] 
below. We will often refer to this weaker version of © as it allows to relax 
some extra conditions imposed on the random environment. 

Condition 2 The distribution of X is non-lattice and belongs without centering 
to the domain of attraction of a stable law with index a € (0, 2] and skewness 
parameter tp, such that 

• \ip\ < 1, ifO < a < I, 

• \ip\<l,ifl<a<2, 

• if) = if a = 1 or a = 2. 

Recall that the condition that the centering constants are zero is redundant, 
if < a < 1, and that a stable law from Condition [2] has the characteristic 
function (see, e.g., Theorem 8.3.2 in [T7]) 

X a ^(t) = exp(-c|t| a (l-#^tan^)), teR, 

where c > is a scaling parameter. This stable law has finite absolute mo- 
ments of all orders r < a. One can check (see, e.g., Section 8.9.2 in [17]) that 
Condition [2] implies ([3]) with parameter 

11 / , ira\ , . 

p=--\ arctan I tp tan — — I G (0,1). (4) 

2 na \ 2 / 

Let 

L. n :=min(5i,...,5„), M n := max (Si, S„) , 

and 

r (n) := min{i £ [0, n] : Si — min (0, L n )} 

be the left -most point at which the ARW attains its minimal value on the time- 
interval [0, n]. It is known (see [3H1 Ch. IV, §20]) that under the Doney-Spitzer 
condition 

n~ 1 r(n) — > r, n — > oo, (5) 

where r is a random variable having a Beta-distribution with parameters 

(1 — p, p) and the symbol -4- stands for the convergence in distribution. 
Let 

OO 

D:=^kT x F(S k =0), 

k=l 
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and 



70 := 0, 7j- +1 := min(n > jj : S n < S 7j ), 

T := 0, := min(n > Tj : S n > S Tj ), j > 0, 

be the strict descending and strict ascending ladder epochs of the ARW. Next 
we introduce two renewal functions 

f 1{0<*} + E£ i P(Sr, <x) if x>0, 

U(x) = ! e-° if x = 0, (6) 

[ if x < 0, 

and 

V(x) = I 1 if x = 0, (7) 

[ if x < 0. 

It is known (see [57], Lemma 1) that (J3|) implies 

E [C/(-X)l {x<0} ] =e-°, 

E[U{x-X)l {x<x} ]=U{x), x>0, 

E[V(x + X)] = V(x), x>0. (8) 



2.3 Refined classification in the subcritical case 

According to (B) the ARW in the subcritical case should have a clear trend 
towards minus infinity. It turns out (see [T], [3T], and [35]) that the behavior of 
a subcritical BPRE to a large extend is determined by the speed of the negative 
drift quantified by the parameter /3, when exists, such that 

E [XeP x ] = 0. (9) 

We will distinguish among three different sub-cases: 

• weakly subcritical, if ([9]) holds with some < /3 < 1, which implies 
E [Xe x ] > 0, 

• intermediately subcritical, if @ holds with j3 = 1, so that E [ATe x ] = 0, 

• strongly subcritical, if E [ATe x ] < 0. 

Notice that due to monotonicity of the function E [Xe' 5X ] in all three sub- 
cases we have E [X] < 0. Clearly, the smaller is the positive value of /3 the closer 
we get to the critical case with E [X] = 0. On the other hand, the subcritical 
case in a constant environment corresponds to X being a negative constant 
implying the strong subcriticality. 

Observe next that E \eP "1 = 7", where 

7 = E[e^] . (10) 
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Let us introduce the following auxiliary measure P with expectation E. For any 
n£N and any measurable, bounded function tp : A™ x Nq 1+1 — > M, the measure 
P is given by 

E[(^Fi, . ..,F n ,Z ,..., Z n )} := [ip{F x , . . . , F n , Z , . . . , Z n )e 0S "] . 

Notice that E[A] = making S a recurrent random walk under P. Clearly, in 
the critical case the measures P and P coincide since (3 = 0. 

Condition 3 The distribution of X with respect to P is non-lattice, has zero 
mean and belongs to the domain of attraction of a stable law with index a S 
(1,2]. 

Under Condition [3] there exists an increasing sequence of positive numbers 
a n regularly varying at infinity 

a n = n Y l a l n (11) 

such that the scaled ARW a~ x S nt P-weakly converges to a strictly stable Levy- 
process L t with parameter a. Here and elsewhere in the expressions like S n t 
the index nt is understood as its integer part. 

3 Annealed approach 

Different trajectories of the ARW represent different scenarios of environmental 
history. The more favorable scenarios translate into higher ARW trajectories. 
If an ARW trajectory stays below zero for a longer period of time, the BPRE 
is doomed to die out. The annealed analysis of the survival of BPRE gives a 
summary picture of the system behavior after averaging over all successful ARW 
trajectories. 

3.1 Asymptotic behavior of the survival probability 

The first result on the asymptotic behavior of the survival probability for the 
critical BPREs has appeared in the seminal paper by Kozlov [42] ■ Kozlov has 
shown for BPREs satisfying Condition [1] that under certain additional assump- 
tions 

cm" 1 / 2 < ¥{Z n > 0) < c 2 n-^ 2 , 

where < c\ < c 2 < oo. Additionally, he proved that in the linear- fractional 
case 

F(Z n > 0) - cn~ 1/2 , c> 0, n -> oo. (12) 

These results were published in 1976. Only in 2000 Geiger and Kersting [3S] 
were able to demonstrate the validity of (| 1 2[) for arbitrary BPREs satisfying 
Condition [TJ To establish this result Geiger and Kersting developed a new 
powerful method of proving conditional limit theorems for BPREs based on a 
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change of measures. This method, along with the idea of splitting trajectories of 
the ARW at the point of global minimum on the interval [0, n] , first suggested 
by Vatutin and Dyakonova [5B], became one of the main tools for obtaining 
conditional limit theorems not only for the critical BPREs but for the subcritical 
BPREs as well. 

3.1.1 Critical case 



The main idea in finding the asymptotic behavior of the survival probability of 
critical BPREs under annealed approach is to show that 

F(Z n >0) ~ 9 P(min(Si,...,S n ) > 0) , n -> oo, (13) 

where > 0. This is done by splitting trajectories of the ARW at the point of 
global minimum on the interval [0, n]. For the annealed approach the trajectory 
splitting method was first used by Dyakonova, Geiger, and Vatutin |30j and 
then further developed by Afanasyev, Geiger, Kersting, and Vatutin [5]. 

Recall that X = In £ [£] is the logarithm of the mean for the conditional on 
the environment offspring number distribution 

Qi = V{ti = i), i = 0,1,2,.... (14) 

In terms of the conditional probabilities Qi the random PGF F can be expressed 
as 

oo 

F{s) = £[st] = ^Q iS \ 

In [5] a higher moment assumption on the environment is given in terms of the 
standardized truncated second moment of the offspring number £: 

C(n) — £ ^ 2l ^ Q >] (15) 
((a) - {s[c])2 ■ (15) 

Condition 4 For some positive e and a 

E [(ln+C(a))^ +e l < oo and E[V{X){ln + ((a)) 1+e } < oo, 

where ln + x := In (max{x, 1}). 

Next we present a number of cases when this assumption is fulfilled. 

1. If the offspring number has a bounded support, i.e. 

P(£ < a*) = 1 (16) 

for some a*, then ((a) = P-a.s. for all a > a* . Obviously, in this case 
Condition |4] is valid. As a particular example we mention here the binary 
splitting reproduction law characterized by 

x 



p(£ = 2 ) = l-p(£ = 0) = e A /2 
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2. In view of relation ©, we have E [V(X)] = V(0) = 1. Therefore, Condi- 
tion [4] is satisfied, if £(a) is a.s. bounded from above for some a £ No. 

Observe that £(2) < 2ry, where 

*"(!) = gfetf-l)] (17) 
(F'(l)) 2 (£[£]) 2 

implying that £(2) is bounded from above if 77 is bounded. This is the case 
when (|14p is either a Poisson distribution with a random mean (when r\ = 1 
a.s.) or a geometric distribution with a random mean (when r\ = 2 a.s.). 

3. The renewal function V(a;) always satisfies V(x) — 0{x) as x — ► oo and 
V(a;) = for x < 0. Therefore, as follows from Holder's inequality, Con- 
dition |4] would be satisfied, provided 

E[(X+) P ] < oo and E [(ln+ C,{a)) q ] < oo 

for some p > 1 and g > max(p _ ,p{p — l) -1 )- 

Under Condition [2] we have p < a -1 due to This observation leads to 
the following weaker version of Condition 01 

Condition 5 For some positive e and a 

E [(ln+C(a)) Q+e ] < oo . 

It was proved in [8] that if ([3]) holds together with Condition|4]or if Condition 
[2] with | tp\ < 1 and Condition [5] are valid, then there exists a positive finite 
number 9 such that, as n — > oo, 

P(Z„>0) - 9F(mm(Si,...,S n ) > 0) - 6 n - {1 - p) l ni (18) 

where ii,/2, ... is a sequence varying slowly at infinity (consult [5] for repre- 
sentations of 9 and the sequence I). Moreover, it was shown in [8] that for a 
trajectory of the ARW to give an essential contribution to the annealed proba- 
bility of survival during the time interval [0,n] this trajectory should pass the 
point of its global minimum over this time interval very soon after its start at 
time 0. 

The result (IT8t gives the asymptotics of the tail distribution for the extinc- 
tion time T = min{fc : Zk — 0} in view of the equality 

¥(Z n > 0) =P(T > n). 

The study of the asymptotic behavior of P (T = n) for the critical BPRE was 
initiated by Vatutin and Dyakonova |53j . They demonstrated that for the linear- 
fractional BPREs under (|3|) , Condition [1] and some additional conditions one 
has 

P(T = n) ~ cn" 3/2 , n-s-oo, (19) 

where c £ (0, oo). Afterwards, Boinghoff, Dyakonova, Kersting, and Vatutin [20] 
have proved (fT9)) for the linear-fractional BPREs under weaker conditions: ([3]), 
Condition [TJ Condition [5] with a — 2 plus the following Condition 
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Condition 6 There exists a constant x € (0, 1/2) such that 

F(x<Qo<l-X, »7>X) = 1. 
Vatutin and Wachtel [51] for the geometric case with 

e -x„ 

F "( s ) = T, — ' n= 1 ' 2 ' 

1 + e — s 

under Condition [2] restricted to a < 2 and |-0| < 1 have shown that similarly to 
(fl"3|) one has 

P(T = n) - 9P(T" = n) as n -> oo, (20) 

where 8 > and 

T - = min{fc > 1 : S k < 0}. 

3.1.2 Subcritical case 

Weakly subcritical case. The best annealed results up to this moment for 
weakly subcritical BPRE were obtained by Afanasyev, Boinghoff, Kersting, and 
Vatutin [11] . It was proved that under Conditions [3] and [Jj (see below) there 
exist numbers < K, k' < oo such that 

P(Z n >0) - KF(mm(Si,...,S n ) > Q) ~ k'— , 

na n 

where 7 is from (fTU)) and a n are from (jlip . Observe that this result is similar to 
that in the critical case (|18p . whereas it is no longer true in the intermediately 
and strongly subcritical cases (sec e.g., 36 ). 

Condition 7 For some positive e and a 

E[(ln+C(a)) Q+£ ] < 00. 

Intermediately subcritical case. Afanasyev, Boinghoff, Kersting, and Vatutin 
[TO] proved that given /3 = 1 and Conditions [3] [7] there are a constant < 6 < 00 
and a sequence l n slowly varying at infinity such that 

P(Z„>0) ~ 07" F(S„ < min(S u ...,S n -i)) rTT^j » ^ °°- 

Under a slightly different assumption the last relation was first proved in |54) . 

Strongly subcritical case. Guivarc'h and Liu in [37] have shown for strongly 
subcritical processes satisfying 

E [£hi + £] < 00 (21) 

that for some < c < 1 

F(Z n > 0) ~ c(E[^])" , 00. (22) 
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This asymptotic formula was originally established by D'Souza and Hambly in 
[48] under an extra moment assumption (also a similar statement was proved 
in [3H] assuming ergodicity of the environment). For the fractional- linear case 
the statement (|22p was obtained in 0]. 

Condition (|21[) is a counterpart of the classical xlogx condition for super- 
critical and subcritical branching processes in constant environment (see |15jV 
It holds, in particular, under the restriction (|16|) . 

3.2 Conditional limit theorems 

Given that a BPRE has survived during a long period of time, what can be said 
about its demographic history? The theorems of this section reveal different 
successful survival strategies adjusted to the variable reproduction strength of 
the population living in a random environment. 

3.2.1 Critical case 

Afanasyev, Geiger, Kersting, and Vatutin [8] proved so far the best annealed 
conditional limit theorem. According to this theorem a surviving critical BPRE 
behaves in a 'supercritical' manner. Supercritical branching processes (whether 
classical or in random environment) grow exponentially fast Z n e~ Sn — > W a.s., 
where W is a non-degenerate random variable (see [TS], @S]). 

For integers < r < n let Xf' n ,t £ [0,1], be the rescaled generation size 
process given by 

Xl' n = Z r+(n _ r)t ■ e -*H-(»-r) t , t £ [0, 1] . (23) 

In [8] it was proved that if (J3)) and Condition [4] (or Condition [2] with \ip\ < 1 
and Condition [5]) are valid, then, 

C(Xt n ' n , t£ [0,1] \Z n >0) ^C(W t ,t£ [0,1]) as n -> oo, (24) 

where the symbol — > stands for the convergence in distribution in the Skorokhod 
topology in space D[0, 1] of cadlag functions on the unit interval or in other 
spaces which we will meet later on. Here r±,r2, ■ ■■ is a sequence of positive 
integers such that r n < n and r„ — > oo. This sequence is introduced to exclude 
the early part of the population history. 

The limiting process {Wt, t £ [0, 1]} is a stochastic process with a.s. constant 
paths, i.e., F (Wt = W for all t £ [0, 1]) = 1 for some positive random variable 
W. We see that the growth of the BPRE is mainly determined by the ARW, 
namely by the sequence (e' s ' l )„>o, with the fine structure of the random envi- 
ronment being summarized by W. Afanasyev |5] in the case r n = established 
the convergence result for t £ (0,1) under a stronger Condition [TJ 
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It was also shown in [5] that given Condition with < 1 and Condition[5] 
there exists a slowly varying sequence l±, 1%, . . . such that 



Z n > 0) 4 £(£+, i G [0,1]) as n oo, (25) 



£(^n o Z„ lnZ„ t , f G [0, 1] 

where L^, t g [0,1], denotes the meander of a strictly stable process with index 
a, which is a strictly stable Levy process conditioned to stay positive on the 
time interval (0, 1] (see [22] and [24]). Convergence (|25l) was earlier established 
by Kozlov [43 under the second moment assumption on the increments of the 
ARW. 

Afanasyev [5] studied the behavior of the process up to the moment of ex- 
tinction and proved that for a critical BPRE satisfying Condition [1] and some 
additional assumptions, asn->oo, 

C (-j= \n{Z nt + 1), t G [0, oo) | Z n > o\ 4 £ (W+ * G [0, oo)) . 

The limiting process W^, t G [0, oo), is a Brownian meander for t < 1 

continued by a " stopped" Brownian motion Wt,t > 1 with Wi = W± killed at 
the moment 

t = inf{t > 1 : W t = 0}. 



Bdinghoff, Dyakonova, Kersting, and Vatutin [20] consider a linear-fractional 
BPRE satisfying Conditions!]] [HI and Condition[5]with a — 2. They have shown 
that 

C(Z n \T = n + l)Ac(Y), rwoo, (26) 

where Y is a positive integer-valued random variable being finite with probabil- 
ity 1, and furthermore, that for any S G (0, 1/2) as n — > oo 

C (z nt e- Snt ,t G [5,1-5] T = n + l) 4£(W t ,i G [(5,1-5]). 

Here the limiting process is a stochastic process with a.s. constant paths, i.e., 
P(Wt = W for all f G [0, 1]) = 1 for some positive random variable W, and 
convergence in distribution holds with respect to the Skorokhod topology in the 
space D[6,l — S] of cadlag functions on the interval [6, 1 — 5]. 

Vatutin and Wachtel [64^ have established for the critical BPREs whose 
increments of the ARWs belong to a-stable law with a < 2 the following 
phenomenon of sudden extinction. Given the process dies at a remote time 
T = n+ 1, the log-size lnZ nt grows roughly as n 1 /" up to moment n and then 
the process instantly dies out, so that \nZ n is of order n 1 /" while Z n+ \ = 0. 
This may be interpreted as the evolution of the process in a favorable envi- 
ronment up to moment n stopped by a sudden extinction of the population at 
moment T = n + 1 due to a "catastrophic" change of the environment. Note 
that this phenomenon is in a sharp contrast with the case of finite variance (see 
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(J2U)) where the distribution of the number of individuals in the population just 
before the moment of its extinction converges to a distribution concentrated on 
the set of positive integers. 

Let us formulate the results of [M] precisely. According to Durrett [H] under 
the Condition [5] there exists a nonnegative random variable A such that for a n 
from (JTTJ) 

lim P{S n < xa n \Si > 0, . . . , S n > 0) = P(A < a:) for all x > 0. 

n— f oo 

It is proved in [64 that if Condition [5] holds restricted to a < 2 and < 1 
together with Condition O then for every x > 

,. P(Z„_i >e a;a ";T = n) E[A- q ;A>.t] 
urn — = B — — : -, 

n^-oo P(T" = n) E [A- Q ] 

where (compare with ([20)) ) 

e:=liminf P(T = n) 
n^oo p(T- = n) 

Moreover, it is found that in the geometric case 

E[A" Q ; A > x] 



UmP(^ 1>e -|T = n)- e[a _ q] . 
3.2.2 Subcritical case 

Weakly subcritical case. Afanasyev, Boinghoff, Kersting, and Vatutin [TTj 
proved for weakly subcritical processes that under Conditions [3] and [7] the con- 
ditional laws C{Z n | Z n > 0), n > 1, converge weakly to some probability dis- 
tribution on the natural numbers. Moreover, it was shown that the sequence 
E[Z„ \Z n > 0] is bounded for any •& < fj, implying convergence to the corre- 
sponding moments of the limit distribution. It was also shown that for Xl' n 
defined in (El 



£(X ( r "'",fe [0, 1] | Z n > 0) 4 £(W t) i G [0, 1]), n -> oo, 

weakly in the Skorokhod space I?[0,1], where ri,r2,... are natural numbers 
such that r n < n/2 and r n ~ > oo. Here Wt = W a.s. for all t E [0, 1] and 
P(0 < W < oo) = 1. Compared to the critical case result now both the 
initial and the most recent part of the population history are excluded. Earlier 
versions of this results can be found in [4] and [36] . 

Thus, we have the following pattern characterising the weakly subcritical 
case: given Z n > the value of Zk are of bounded order for k close to and 
close to n. For 1 <C k -C n the demographic dynamics of Z^ follows the value 
of the quenched conditional mean £ [Zk] = e Sk in a completely deterministic 
manner up to a random factor W > 0, resembling the behavior of supercritical 
branching processes. 
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Intermediately subcritical case. For intermediately subcritical BPREs 
under the annealed approach Afanasyev, Boinghoff, Kersting, and Vatutin [10] 
proved that under Conditions [3] and [7J the unsealed population size Z n condi- 
tioned on {Z n > 0} converges in distribution asn-> oo. Under slightly different 
assumption this result was first obtained by Vatutin in |54) . 

Denote by L t the strictly stable process L t with parameter a conditioned to 
have its minimum at time 1. Let e-j,, eg, . . . denote the excursion intervals of L t 
between consecutive local minima and put j (i) = i for t € ej. In 10 it was also 
proved that Conditions [3] and [7] imply for <t\ < t 2 < ... < tk < 1, as n — > oo, 

£ f Znh Z ntk 

\exp (Snti - nainfe< ntl S k ) ' ' exp (S n t k - min fe < Iltfc Sk) 

where Vi,p2>--- are i.i.d. copes of some strictly positive random variable V. 
This means that if ti and tk belong to one and the same excursion, then j (ti) 
and j (tk) coincide, i.e. VjUA = Vju k ) w hh probability 1. Notice that here 
j(ti) _• j(*2) < ■ • ■ < j(t n )- For the linear fractional case this result was earlier 
obtained by Afanasyev in [7J. 

Recall that r (nt) is the time of the first minimum of the ARW up to time 
[nt] . Complementing the previous results it was shown in [10] that at the times of 
consecutive ARW minima the population sizes have discrete limit distributions 

£( (Z T(ntl)) . . .,Z T(ntk) ) | Z n > ) 4 C (Y, (tl ), . . .,Yf(t*,)) , 
where Y\, Y%, . . . are i.i.d. copies of a random variable Y taking values in N. 

Strongly subcritical case. Geiger, Kersting, and Vatutin [36] studied the 
strongly subcritical BPRE under condition ([2T]l and proved that 

lim P (Z n = z I Z n > 0) = r z , z £ N, 

n— >oo 

where 

oo oo 

r z = 1 and m r := zr z < oo. 

2=1 2=1 

Afanasyev, Geiger, Kersting, and Vatutin [9] have introduced the following 
condition in the strongly subcritical case. 

Condition 8 Suppose that (see (|17l) ) 

E [e x ln+ -q] < oo. 

Since in the strongly subcritical case E [e x ] < oo, Condition [S] holds in par- 
ticular if the random offspring distribution (| 14[) has uniformly bounded support. 
It also holds if ([14]) is a Poisson distribution with random mean, so that rj = 1 
a.s., or if ([M]) is a geometric distribution on Nq where rj = 2 a.s. 



Z n > 
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It was shown in [pj that with 

=: i n ,o < Vi < i n ,2 < ■■■ < i n ,k < Vfc+i : = n 
the following weak convergence holds 

^( (^i)o<i<m ' i^in.i+o) 0<i<m ' "•' (^n,*+i)o<j<m ' ^ n -^0<j<m I ^" > ^) 
^ («) <,< ro («) <i<™ f*®^ {(Xi)0<i<m ), 

for every k,m £ No as mino<z<fe («n,z+i — *n,i) — >• co and n — > oo. Here, 
£„ ((^j)j> ) denotes the law of the Markov chain (Yj)j >Q with initial dis- 
tribution v and transitional probabilities 



p _ *p(ft + -+ey = *) N 



where £1, . . . ,£ a are independent copies of the offspring number £. The station- 
ary distribution of {Yj)- >0 is 

r z = , z G N. 

The process (Xj)j>o is the time- reversed Markov chain with transition proba- 
bilities ^ 

'2/ 



4 Multi-type BPREs 

The model of multi-type BPRE was first considered by Athreya and Karlin [T3] 
and subsequently investigated by Weissener [ST], Kaplan [3D], Tanny [?pj ( see 
also for other references). Next we describe the BPREs with p types of 
particles Z„ = {Z n> \, Z n ^ p ) and time n = 0, 1, . . . using standard notation for 
the p-dimcnsional vectors: 

• the unit vector has j-th component equals 1 and all others equal zero, 
3 = >P) 

• all zero and all one vectors = (0, . . . , 0), 1 = (1, . . . , 1), 

• for x = (xi,...,x p ) and y = (yi, ...,y p ) set 

p p 

ll x ll = X^l^l' ( x >y) = x: 

i=l »=1 



=n-r 

! = 1 
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A particle of type i in a multi-type branching process can produce z\ par- 
ticles of type 1, Z2 particles of type 2, . . ., and z p particles of type p with a 
probability Qz \ z = (21, ... , z p ). Under the random environment assumption 
the probabilities Q% are random variables making random the corresponding 
PGFs 

OO OO 

F«( Sl ,..., Sp )=]T...^QWs z , i = l,...,p. (27) 

21—0 z p — 

The random PGFs for the consecutive generations (F„ , F^) n >i represent- 
ing the environmental history determine the conditional mean offspring numbers 
as random p x p matrices M„ with elements dF^f 1 (l)/dsj. These are IID ma- 
trices having the same distribution as 

dF^(l) 

M = (M«)? i=1 , Mij = —^-. (28) 

OSj 

We will assume that elements My are all positive meaning that each of p 
types of particles can produce any other type of particles in the next generation. 
We denote by R and R n the Perron roots, i.e. the maximal (in absolute value) 
eigenvalues, for the matrices M and M„. There exist left and right eigenvectors 
corresponding to the Perron eigenvalue 

V=(Vi,...,V p ), VM = i?V, 
U = {U X ,...,U P ), MU' = i?U', 

both with positive components and unique up to scaling. Although vectors V 
and U are random in general, throughout this section we assume that one of 
the following basic conditions is valid. 

Condition 9 There is a non-random strictly positive vector v such that ||v|| = 
1 and 

P(V = v) = 1. 

Condition 10 There is a non-random strictly positive vector u such that ||u|| 
1 and 

P(U = u) = 1. 

The important gain of this restriction is that we get either 
vMi . . . M„ = R! ■ ■ ■ i?„v 



or 

Mi . . . M n u' = Ri ■ ■ ■ R n u'. 

In both cases putting X := lni?, X n := lni?„, n > 1, we can again introduce an 
ARW by S n — X\ + ... + X n . We will refer to the conditions on S n stated for 
the single type case even in the multi-type setting (under Conditions l9l or [TO]) . 
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In particular, the single type classification for BPREs extends straightforwardly 
to this multi-type case. For example under Condition [9] as n — > oo 

p 

V ViV{Z n I Z = e,) < min |vMi • • • M fe | = min i?x • • • R k 

^— ' Kk<n Kk<n 

i=l - - - - 

= exp{ min /Sfcl- — ^ 

l<fc<n 

almost surely. This implies that in the critical and subcritical cases 

V(Z n f I Z = ei ) -> 

almost surely for any i = 1, . . . ,p. The same conclusion follows from Condi- 
tion rTU] as well. 

4.1 The survival probability of the multi-type BPREs 

Dyakonova and Vatutin have managed to extend some of the single type an- 
nealed results to the multi-type BPREs in the critical and subcritical cases. 
Some of the results mentioned in this section bring a stronger version of the 
known single type results. 

4.1.1 Critical multi-type case 

Dyakonova [28] considered a multi-type critical BPRE satisfying ([3]), Condi- 
tion [SJ and Condition @] where ((a) is now defined as 

C(a):=iT 2 J2 I> E Q * Z ^ ^ 

zSNg i=l j,fc=l 

with being the set of vectors z with non-negative integer- valued components 
such that at least one component is larger or equal a. Here R is the (random) 
Perron root of the random matrix M related to through (|2~T|) and (|2"5)) . 
Additionally the following restriction (due to [H|) is imposed. 

Condition 11 There exists a number < d < 1 such that 

d< — ^ <<T\1 < h,i2,h,32 <P- 

Under these conditions it was proved in [25] that (compare with l[T8|) ) 

P(Z„ + | Z = ei ) ~ an-^-Phn, i = l,...,p, (30) 

where Ci are positive constants, p is the constant from the Spitzer condition 
applied to S n — lni?i + ... + mi?„, and l n is a sequence slowly varying at 
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infinity. Earlier the same asymptotics (f3U)) was established in [26 in the linear 
fractional case under Condition Q] with \ip\ < 1 and Condition [TUl 

In [35] it was proved also a quenched type result claiming the almost sure 
convergence 

V(Z n ? I Z = e !: ) 

> U , z = l,...,p, (31) 



E?=i«i^(Z»^O|Z o =e J 

where U* = ([/-j", t/*) is a random vector such that 

(v,U*) = 1, d< [/; < l/min^,...,^), i = l,..., p. 

Furthermore, it was shown in |28j that the above mentioned results remain true 
for an alternative set of conditions involving the counterparts of the single type 
Conditions [2] and [5] 

Vatutin and Dyakonova [52] studied a multi-type critical BPRE whose ARW 
satisfies the Spitzer-Doney condition and the mean matrices of the reproduction 
laws have a common positive right eigenvector (see Condition HUl) . Also in [62] 
it is assumed that the following assumption holds. 

Condition 12 By this condition we exclude a possibility for the defining PGFs 
to take the linear form 

F®(a) = Q®+Q®8 1 + ... + Q®8 p . 

Another condition required in |62] concerns the moments for the vectors 
(Sti) • • • i £,ip) °f the offspring numbers having quenched distributions {Q { : } } and 
quenched means (Ma, . . . , M ip ). For a given /3 > set 

Aij {(3) = £ |& - Mij f , A = max Ay {(i) 

ho 

and denote C,p := e~@ x Ap. 

Condition 13 There exist (3 G (1, 2] and e > smc/i i/iai 

P(A,g < oo) = 1, (32) 

and 

E[ln+ C^] 1/p+£ < oo, E[V(X)(ln+ C^) 1+£ ] < oo. 

In [62] the asymptotic result (|30]l was established under ([3]) and Conditions 
ITOl [T2l and HU The results of [62] generalize and extend not only the above 
mentioned results in [35] but also the corresponding single type statements from 
[5J . Indeed, in view of condition (|32|) the asymptotic relation (|30l) may be valid 
for p = 1 if, for instance, 

¥(£ [ei {i > a} ] =oo) = l 
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for any a > making in this case the conditions of the last result weaker than 
those required in [8] for (fT8|) . 

As in the single type case, denote by T the extinction moment for the process 
Z„ . It was shown by Dyakonova [27] that in the multi-type linear fractional case 
under Conditions IT1 l9l HU1 ITT1 and some extra restrictions there exists a constant 
< c < oo such that (compare with (|19|l ) for any i = 1, . . . ,p 

P(T = n\Z a = e 4 ) - c Ul n~ 3/2 , n ->■ oo, 

where u = (tti, u p ) is defined in Condition 1101 

4.1.2 Subcritical multi-type case 

Dyakonova [29] investigated the asymptotics of the survival probability of multi- 
type subcritical BPREs. She has proved, that if Conditions [9l and [TT1 are valid 
and E [X] < 0, then ([31]) holds P-a.s., as n — > oo. Furthermore, in [29] it was 
shown that in the strongly subcritical case as n — > oo (compare with ([22])) 

P(Z„ ^ | Z = e 4 ) - q(E [i?])", a > 0, i = 1, . . .,p 

under Conditions [3 HU and [8] with rj = C(0) defined by (p9|) . 

4.2 Functional limit theorem in the critical case 

Vatutin and Dyakonova 62| have assumed that assumption ^ is valid together 
with Conditions [TU1 [T2l and Condition [T51 with /3 = 2. They considered a family 
of the processes 

W r , n (t) :=e- s H-c»-r )*(Z r+(n _ r)t ,u), te [0,1] 

for any given pair of integers < r < n. Letting ri,r2, . . . be a sequence of 
positive integers such that r n < n and r n — > oo, it was proved that as n — > oo, 

C(W rn , n (t),te [0,l]|Z„^O;Z o = z)4£(VF z (0,te [0,1]), (33) 

where the limiting process is a stochastic process with a.s. constant paths, i.e., 
V(W z (t) = W z for all t G [0, 1]) = 1 with P(0 < W z < oo) = 1. Convergence 
([33]) is a generalization of (f24|) in two directions. First, it is proved for the 
multitype case, and second, even for the single-type case the conditions under 
which ([33)) is established in are weaker than those used in [5]. 

5 Quenched results for the critical BPREs 

Vatutin and Dyakonova [56 and 57 applied the quenched approach to the crit- 
ical BPREs meeting the Doney-Spitzer condition ([3]). Their results generalize 
those in [55] and [58] in which much stronger conditions on the characteristics 
of the branching processes are imposed. 
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5.1 Survival probability conditioned on the environment 

The following condition involving the renewal functions © is dual to Condition 
2] dealing with the renewal function (J7J . 

Condition 14 There are numbers e > and a £ No such that 

E[(ln + C(a))^ +E ] < oo and E[[/(-X)(ln+ C(a)) 1+ 1 < oo, 
where £(a) is /rom L?5j) . 

It was proved in [57] that under ([3]) and Conditions [4j [14] 

e - 8 rwp (Z n > 0) 4 C, n -> oo, (34) 

where the random variable Q € [0, 1] is positive with probability 1. According 
to l)34p the asymptotic behavior of the survival probability is mainly determined 
by the minimal value of the ARW over the time interval [0,n]. 

5.2 Convergence of one-dimensional distributions 

Denote 

Y ■= — 

n ' £ [Z„ | Z n > 0] 

and 

X (n) (da;) := V (? n e |Z„ > 0J . 

Theorem 1 and Lemma 7 from [57] yield that under ([3]) and Conditions U] 
1141 there exists a random measure Ai (being proper and nondegenerate with 
probability 1) such that, as n — > oo, 

/•oo 

/ e- Xx M (n) {dx)A<$>{X) 
Jo 

POO 

/ e- Xx xM (n) (dx) 4 *(A) 
Jo 

In the linear-fractional case the limiting Laplace transform is deterministic 
$ (A) = (1 + A) -1 and represents an exponential distribution with parameter 
1. Notice that in this case ^ (A) = (1 + A)~ 2 corresponds to the sum of two 
independent Exp(l) random variables. Thus in the linear-fractional case we get 
an analogue of the corresponding Yaglom-type limit theorem for the ordinary 
critical Galton- Watson processes. 

For the case when the distribution of X is absolutely continuous convergences 
(|34|) and ([35]) were proved in [58] . 



e- Xx M(dx), 
e~ Xx xM{dx) 



d\ 



(35) 



$ (A) . (36) 
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5.3 Convergence of finite-dimensional distributions 

For an integer b > 2 and tuples t — (to, ti, . . . , 4), A = (Ai, . . . , A;,), where 
= to < *i < " • • < *6 = 1 an d Aj > 0, i = 1, . . . , b, set 

$ (n) (t,A) :=£ 

and for a positive integer d < b and a vector r = (ro, ri, r"2, • • • , Td) with integer- 
valued coordinates = ro < T\ < T2 < ■ ■ ■ < — b introduce the event 

U(i;r,n) := {r(raii) = • • ■ = r(nt ri ) < r(nt ri+ i) = • • ■ 

= r(nt r2 ) < • ■ ■ < T{nt rd _ 1+1 ) = ■■■ = r(nt r j|. 

It was proved in [57] that if X meets Conditions [2] [U [141 then, as n — » oo, 

d-i ( Ti \ ( b 

i=l \j'=r<-i+l / \j=r d -i+l 

where ^^A), i = l,...,d— 1 and $d(A) are independent random functions 
distributed as ^(A) in (1551) and $(A) in (|55|) respectively. 

Loosely speaking, this result shows that the trajectory (Z , ...,Z n ) is parti- 
tioned into independent pieces generated by strict descending ladder moments 
of the ARW. The population size at the moment j is proportional to e S 3~ s tU) 
while the maximal size of the population on the interval [0, n] is proportional to 
e max < J <„(s J -S T(i) ) (with a random positive factor < 1). 

Let L t be an a-stable Levy process. For fixed = to < ti < • • • < t& = 1 
put Up := inf to < u <4 p L u , < p < b. Let D be the random number of different 
elements in the tuple u>q, w\, . . . , u>b, and let R be the random Z?-dimensional 
vector with components 

Ro = 0, R i+ i = max {k > R4 + I: u Ri+ i = co k } 

so that Rd — b. In [57] it was established that if X satisfies Conditions [2] [4] 
[LA then as n — > 00 

$w(t,A)4nV e e a 

where (\&i(A), $i(A), ^2(A), $a(A), . . .) are as before and independent of the 
Levy process L t defining the vector R. 

Relations and (|55|) imply that if t G (0, 1] is fixed then, given Z n > 0, 
the distribution of the random variable Z nt e ST{nt) ~ Snt converges in the specified 
sense to a proper distribution with no atom at zero. This means, roughly 
speaking, that if the process survives up to moment n, then an earlier population 




Z n > 
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size Z nt is proportional to e Snt ~ Srf - nt '> . Thus, contrary to the conditional limit 
theorems for the classical critical or supercritical branching processes (in which 
the scaling functions for the population size increase with time either linearly or 
exponentially), the scaling function for the population size in the critical BPRE 
is subject to large random oscillations. 

This indicates that the corresponding population passes through a number 
of bottlenecks at the moments around the consecutive points of minima of the 
ARW. Vatutin and Dyakonova in [55J, [53], [57], [55], and [5j|] have investigated 
this phenomena in detail under the quenched approach. They shown that the 
distribution of the number of individuals in the process at the moments close 
to r(nt), t £ (0, 1] conditioned on survival up to time n converges to a discrete 
distribution. Thus, in contrast to the fixed points of the form nt where the size 
of the population is exponentially large (see [57] ), the size of the population at 
the (random) point of global minimum of the ARW becomes drastically small 
but then it grows again exponentially. This reminds the typical demographics 
of real biological populations which during their evolution have "favorable pe- 
riods" (rapid growth of the population size) and " unfavorable periods" (quick 
extinction when only a few representatives of the population survive) followed by 
another period of rapid growth. Note a similarity between this picture and the 
intermediately subcritical case under the annealed approach (see Section l3.2.2[) . 



5.4 Discrete limit distributions 



Consider the population size near the time of the global minimum of the ARW. 
Under (J3J and Conditions H] Q3] it was proved in [59] that for any m £ Z and 

te(o,i) 



CIS 



£ 

Z n > f) 



Z n > 

r(rz.) > nt 



where (p m (s) and ^(s) can be obtained precisely. Observe, that the second 
limiting PGF does not depend on a particular value of t £ (0, 1). 

Next we give another two convergence results from [59) . The first of them 
characterizes the conditional distribution of the number of individuals in the 
population at two sequential moments located in a vicinity of the moment when 
the global minimum of the associated random walk is attained, while the second 
convergence describes the same distribution in a vicinity of the moment r(nt) 
given r(n) > nt. So, in [5J5] proved that under ^ and Conditions HI [TH for any 
m £ Z and any t £ (0, 1) , s\, s 2 £ [0, 1), as n — > oo, 



^Z(T(n)+m) Z{T{n)+m+l) 



CIS 



^Z{r{nt)+m) Z(r{nt)+m+l) 



Z(n) > ->• (p m {si,s 2 ), 
Z(n) >0 T{n)>nt\ A C[y* m (s l ,s 2 )) 



where ip m (s±,S2) and <^(si,S2) can be obtained precisely. 
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5.5 Conditioning on the precise time of extinction 

Critical BPREs conditioned on extinction at a given moment were investigated 
by Vatutin and Kyprianou [53J using the quenched approach. Let 



t(L n) := min < k € \l,n] : Sk = min S v 

{ l<p<n 

be the left- most point on [I, n) at which the minimal value of the ARW on 
the interval [l,n] is attained. Assuming {T = n} the authors of [53J proved 
conditional limit theorems describing the asymptotic behavior, as n — > oo, of 
the distribution of the population sizes at moments nt, t £ (0, 1) and at moments 
close to r(nt). It turned out that if X belongs to the domain of attraction of a 
stable law with parameter a £ (0,2], then (contrary to the annealed approach 
[64]) under the quenched approach the phenomenon of sudden extinction does 
not occur. 
Let 

7>(Z„>0|Z = 1) 



T(Z n >0\Z m = l 



-P(Z n =0\Z m = l)b 



in i 



where b m := X^lo Vj+i e /2, and Vj are from (fi~7j) . It is shown in [63] that 
in the linear- fractional case given assumption (J3J and Conditions 2] Q3] are valid 
the following is true. For any m e Z, t e (0, 1), s e (0, 1], A £ (0, oo), as n — > oo 



C (£ [a^(»*)+»|T = n] |r(n) > nt) Ac(s(J- 




(37) 



1 - . 



C (£ [s z ^,n )+m \ T = n ] |r(n) < nt ) 4 £ ^ S ^— ^- J J , (38) 
where 8 m € (0, 1) , 9 m G (0, 1) with probability 1, and furthermore 



£ 



(1+A) 2 



(39) 



According to [53] on the set r(n) > the random variable An^ne^'"*) -5 " 4 
converges in distribution, as n — > oo, to a random variable being finite and 
positive with probability 1. Thus, for such moments nt the normalization in 
(l39|) is essentially specified by the past behavior of the ARW. On the other hand, 
given r(n) < nt the random variable A n t, n e S ' r{nt,n) ~ Snt converges in distribution, 
as n — > oo, to a random variable being finite and positive with probability 1. 
Thus, for such moments nt the scaling in (|39p is essentially specified by the 
future behavior of the ARW. 

This fact allows us to give the following non-rigorous interpretation of the 
mentioned results. If the process dies out at a distant moment T — n, then it 
happens not as a unique catastrophic event. Before the extinction moment the 
evolution of the process consists of a number of "bad" periods characterized by 
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small population sizes. According to (|57|) and (j3"5j) . such periods are located in 
the vicinities of random points T{nt)li T i n )> n t} and r(nt, n)ls T i n \ <nt \. On the 
other hand, at nonrandom points nt, t £ (0, 1), the size of the population by (|39p 
is big. Hence lnZ nt grows like S nt — SV(nt) if r(n) > nt and like S^t — SVfnt.n) 
if r(n) < nt. Thus, the process {Z t T,0 < t < 1} dies by passing through a 
number of bottlenecks and favorable periods. 



6 Reduced BPREs 

The next stage in studying branching processes is to investigate the structure 
of their genealogical trees. For < k < m let Zk. m be the number of particles 
at moment k in the process {Z n , n > 0} each of which has a nonempty offspring 
at moment m. The tuple {Zf. jm ,0 < k < m < oo} is called the reduced branch- 
ing process which gives the number of branches in the genealogical tree of the 
population modeled by the branching process {Z n ,n > 0} . 



6.1 Annealed approach 



The first results for reduced BPREs with linear-fractional PGFs and under 
the annealed approach were established by Borovkov and Vatutin [19] for the 
critical case and by Fleischmann and Vatutin [34] for all three types of subcritical 
BPREs. For the critical case Vatitin [52] proved under Condition Q] and some 
additional conditions that 



C 



1 



; In Z„ 



t e [o, l] z n > o) 4 c ( inf w+, t e [o, 1]) , 



(40) 



where W t , t G [0, co), is the Brownian meander and the convergence in dis- 
tribution holds in Skorokhod topology in space D[0, 1]. This convergence was 
established in Q15] for the case of linear-fractional BPREs under stronger mo- 
ment assumptions. 

Let B n = max{m < n : Z m ^ n = 1}. The difference d n = n — B n is called the 
time to the most recent common ancestor for all individuals existing at time n 
(MRCA„). For ordinary Galton- Watson processes conditioned on survival, the 
classical results by Zubkov [55] state that in the critical case d n is asymptotically 
uniformly distributed over [0, n] while in the subcritical case d n is asymptotically 
finite. We have a quite different situation in the random environment setting. 

For the annealed approach convergence (|40|) states, roughly speaking, that 
the number of individuals Z„ tj „ grows as exp{y/n mi t < u <i Recalling that 

P(inf t <„<i W+ > 0) = 1 for any t £ (0, 1] , we conclude that in the annealed 
setting conditionally on survival the MRCA n has lived at early time of order 
o(n). In fact, as shown by Borovkov and Vatutin [TH] for the linear- fractional 
case, the MRCA n with positive probability is the initial individual! 

Studying the subcritical BPREs with geometric offspring distributions Fleis- 
chmann and Vatutin [34j have found that for the intermediately and strongly 
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subcritical processes the MRCA„ is located not too far from the point of obser- 
vation n, while for the weakly subcritical case it may be located either nearby of 
the point of observation or at the very beginning of the evolution of the process. 

6.2 Quenched results for the critical reduced BPREs 

This section contains a summary of results established by Vatutin and Dyakonova 
in [60] and [61] for the critical reduced BPREs under the quenched approach. 
The results formulated in this section are proved under ^ and Conditions @] 
[Ml if not explicitely stated otherwise. 

6.2.1 Time to the MRCA and afterwards 

Vatutin and Dyakonova [60] proved that in the critical case the MRCA„ is 
located "not too far" from the moment of the global minimum r(n) of the 
ARW on [0, n] in that for any m E Z 

oo 

V(B n = r(n) +m) -+r m , ^ r ™ = 

m— — oo 

and with probability 1 

lim lim V (\B n - r(n)\ < m) = 1. 

m— >oo n— >oo 

Recalling ([5]) we clearly see the difference with the ordinary Galton- Watson 
critical branching processes where the time to the most recent common ancestor, 
scaled by n, is uniformly distributed on [0,1]. 

Also in [60) it was proven that the finite-dimensional quenched distributions 
of the process {Z T ^ +mn , m£Z} conditioned on Z n > weakly converge to the 
finite-dimensional distributions of a Galton- Watson branching process evolving 
in an inhomogeneous random environment. The limit process starts at — oo 
by a single individual. The founder of the population dies at a moment m € 
(— oo, +oo) with probability r m producing at least one offspring in accordance 
with the PGF a m (s) defined as the weak limit 

£ [s z ^ )+m+1 - n \z T(n)+m ^ n = l]A a m (s),se [0,1), 

as n — > oo. The next generation particles reproduce independently according to 
the PGF o m +i (s). And so on. Observe that this description is in a sharp contrast 
with the respective limit process for the ordinary critical reduced Galton- Watson 
processes obtained in [33] . 

The properties of reduced processes far to the right of r(n) were studied in 
[60] as well. Using the scaling function 

, n { 1 £[Z k , n \Z n >0]' ~ ~ ' 



2G 



it was shown that as n — > oo and m — > oo with r(n) + m < n 



£ 



-AZ T( „ ) + m ,„/3„(r(n)+m) 



z n > 



3(A), 



where $(A) is the same random Laplace transform as in (|35|). Moreover, if 
< mi < TO2 < ... < trip < n — T(n), then as mi — > oo and n — > oo 



£ 



IT 

.i=l 



-\iZ T 



(„)+m 4 ,n^(r(n)+m 4 ) |^ > q 



4 * a, 



6.2.2 Reduced processes at nonrandom times 

Properties of the reduced process Z nt , n for nonrandom t £ (0, 1) were also 
described in EH1 ■ It was shown that 



£ 



z„ > 



$(A)l{ T < t } + e A l{r>t}, n^oo 



where r is a random point having a Beta(l — p, p) distribution as in ([5]) and 
being independent from $(A). Furthermore, for any Ai > 0, . . . , Afc > and 
< ti < ... < tk < 1, as n — > oo, 



exp < - ^ KZnU.nPn (nU) \ \ Z n > 



(=1 



4* 5> 



where i r = min{i : r < tj }. 

The corresponding limiting process {W t , < t < 1} follows a simple pattern 
(conditioned on a given environmental development). Until time r (correspond- 
ing to the time to the MRCA) we have Wt = 1 and after time r we have Wt = Y, 
where the random value Y has $(A) as its Laplace transform. Recall that in 
the linear- fractional case Y is exponentially distributed with parameter 1. 

Concerning the scaling function f3 n (nt) it is known from |60j that under 
Conditions HHElfor any t € (0,1) 

C (e s r^t,n)-s^ n) p n (nt)\ T (n) < nt ) 4 C {(*) , n -> oo, 

where C* is a proper random variable. It follows that if the global minimum 
time r(n) of the ARW occurs prior to nt then Z nt: n is of order e t(™*"0 — 
(up to a random multiplier separated from zero and infinity). 



6.2.3 Reduced BPREs as random fields 

Vatutin and Dyakonova [61 obtained quenched limit theorems for the random 
fields { 

}o<ti<t 2 <i as Z n > and n — > oo. The answers essentially 
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depend on the position of r(n) relative to the interval [nt\, nt 2 \. as n — > oo 
C (£[e- xz ^,*Mtut 2 )\ Zn > o] I r(n) < nii ^ 4> $ (A) , (41) 

^^[e-^"*!'"^^^ 1 '* 3 ^,, > 0] r(n) > nt 2 ,r(nt 2 ) < nil) 4*(A), (42) 

£ ^[ s z »*i.™*2 |z„ > 0] < r(n) < nf 2 ) 4 s, (43) 

C (£[s z »*i.»*2 |Z„ > 0] r(n) > nt 2 ,r(nt 2 ) > nt^j 4- s, (44) 
where $ (A) is from ((35]). * (A) is from ([36]). and 

o /. . s ._ 1 1 - /o,Kt 2 (0) 

Pn[1 ' 2) ■ £[Z ntlM2 \Z nt2 > 0] e s -i (l-/„ tl ,nt a (0))' 

According to Theorem 1 in 56 

/3„(ii, i 2 )e S ^" t i'" t ^ SV < n *i ) 4 9, rwoo, 

where the random variable O, is positive and finite with probability 1. Observe 
that in the cases (|43|) and (|44| we have S T ! ntltn t 2 ) < S T r ntl ), Results (|43|) 
and (|44p say that in these cases all individuals existing at moment nt 2 are 
descendants from a single individual existing at moment nt\. Relations (|41l) 
and (f4"2"]) describe the cases when the corresponding number of ancestors is large 
with Z ntunt2 x e s -<"*i,-t 2 )-S T( „ tl) ^ 
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